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ABSTRACT 

We combine constraints on galaxy formation histories with planet formation models, yielding 
the Earth-like and giant planet formation histories of the Milky Way and the Universe as 
a whole. In the Hubble Volume (10^^ Mpc^), we expect there to be ~ 10^** Earth-like and 
^ 10^*^ giant planets; our own galaxy is expected to host ~ 10® and ~ 10**^ Earth-like and giant 
planets, respectively. Proposed metallicity thresholds for planet formation do not significantly 
affect these numbers. However, the metallicity dependence for giant planets results in later 
typical formation times and larger host galaxies than for Earth-like planets. The Solar System 
formed at the median age for existing giant planets in the Milky Way, and consistent with past 
estimates, formed after 80% of Earth-like planets. However, if existing gas within virialised 
dark matter haloes continues to collapse and form stars and planets, the Universe will form 
over 10 times more planets than currently exist. We show that this would imply at least a 
92% chance that we are not the only civilisation the Universe will ever have, independent of 
arguments involving the Drake Equation. 

Key words: planets and satellites: terrestrial planets, planets and satellites: gaseous planets, 
galaxies: formation 


1 INTRODUCTION 

Early estimates of the planet formation history of the Universe 
(Livio 1999; Lineweaver 2001) suggested that the Earth formed 
after 75-80% of other similar planets, even when considering po¬ 
tential galactic habitable zones (Lineweaver et al. 2004). Since that 
time, thousands of exoplanets have been found, aided by the Kepler 
mission (Lissauer et al. 2014). Many advances have been made in 
the past decade, especially in our understanding of how planet for¬ 
mation depends on the mass and metallicity of the host star (Eischer 
& Valenti 2005; Buchhave et al. 2012; Wang & Eischer 2013; Lis¬ 
sauer et al. 2014; Buchhave et al. 2014; Gonzalez 2014; Reffert 
et al. 2015). Concurrently, constraints on galaxies’ star formation 
and metallicity histories have been improving rapidly (Maiolino 
et al. 2008; Mannucci et al. 2010; Moustakas et al. 2011; Behroozi 
et al. 2013e; Peeples & Somerville 2013; Munoz & Peeples 2015). 

In this paper, we combine recent planet frequency models 
(Lissauer et al. 2014; Buchhave et al. 2014) with reconstructed 
galaxy formation histories (Maiolino et al. 2008; Behroozi et al. 
2013e) to update constraints on the planet formation history of the 
Milky Way and the Universe as a whole, both for Earth-like plan¬ 
ets and for giant planets. We adopt a flat, ACDM cosmology with 
SIm = 0.27, = 0.73, h = 0.7, tis = 0.95, and ctg = 0.82, simi¬ 

lar to recent WMAP9 constraints (Hinshaw et al. 2013); the initial 
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mass function (IMF) is assumed to follow Chabrier (2003) from 0.1 
to 100 M 0 . 


2 METHODOLOGY 


While planet formation can depend on host star mass, new stars’ 
masses are drawn from a nearly universal distribution (Chabrier 
2003). When averaged over an entire galaxy, the planet formation 
rate (PER) is then proportional to the galaxy’s star formation rate, 
modified by the PFR’s metallicity dependence. Using a power-law 
parametrisation for this metallicity dependence (see, e.g., Fischer 
& Valenti 2005; Gonzalez 2014), we model the planet formation 
rate of a galaxy as a function of its stellar mass (M* (r)) and cosmic 
time (f): 


PFR(M*,t) 


(Z(M„t)/ZQ)« 

^ / \ 

(m*) 


SFR(M*,r), 


( 1 ) 


where n is the mean number of planets formed per star, a is the 
power-law dependence of planet incidence on metallicity, Z{Mt,t) 
is the galaxy’s mean gas-phase metallicity, (ra*) = 0.67 M 0 is the 
mean mass of a newly-formed star (Chabrier 2003), and SFR(Af* ,t) 
is the galaxy’s star formation rate in M 0 yr^*. Additional factors 
(e.g., stellar initial mass functions and densities) influencing the 
PER are discussed in Appendix A. 

For giant planets (R > 6 R 0 ; including, e.g., Jupiter and 
Saturn), the metallicity (specifically, [Fe/H]) dependence is long- 
established (Fischer & Valenti 2005); recent estimates suggest 
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Figure 1. Top-left panel: formation rate (in planets/yr) for Earth-like planets as a function of galaxy stellar mass and cosmic time. The dashed line indicates 
the median expected growth history of the Milky Way (Behroozi et al. 2013e). The dot-dashed line indicates [Fe/H]=—1.5, which has been suggested (Johnson 
& Li 2012) as the threshold metallicity for planet formation. Grey shaded areas indicate where galaxies are not expected to exist in the observable Universe. 
Top-right panel: same, for giant planets. Bottom-left panel: Earth-like planet formation rate multiplied by galaxy number density as a function of stellar 
mass and cosmic time, i.e., the volume density of planet formation (in planets/yr/comoving Mpc^/dex). Contours indicate where 50% and 90% of all planet 
formation has taken place. The 0 symbol indicates the Milky Way’s stellar mass and age at the formation of the Solar System. Bottom-right panel: same, for 
giant planets. 


hq ~ 0.022 and ac ~ 3.0 (Gonzalez 2014), albeit with signifi¬ 
cant systematic uncertainties. To define Earth-like (i.e., “habitable 
zone”) planets, we adopt the same definition as Lissauer et al. 
(2014), requiring that planets with an Earth-like atmosphere could 
support stable surface reservoirs of liquid water. Effectively, this 
includes all objects whose radii and orbital periods are within a 
factor of e of those of the Earth (Lissauer et al. 2014). In the So¬ 
lar System, this would include Mars and Venus, but exclude, e.g.. 
Mercury and the Moon. The metallicity dependence for Earth-like 
planets is believed to be smaller than for giant planets (Buchhave 
et al. 2012; Campante et al. 2015), with recent estimates suggesting 
~ 0 to 0.7 (Wang & Fischer 2013; Lissauer et al. 2014). This 
range of «£ has only a small impact on our results, so we conser¬ 
vatively take «£ = 0 (i.e., no metallicity dependence) for Earth¬ 
like planets. However, Johnson & Li (2012) suggest a theoretical 
minimum metallicity threshold for Earth-like planet formation of 
[Fe/H]~ -L5+log[o (au) (with r the orbital radius), so we mark 
a fiducial threshold of [Fe/H]= —1.5 in all relevant plots. For he, 
Kepler has provided the largest statistical samples (Catanzarite & 
Shao 2011); Lissauer et al. (2014) suggest an incidence of ~ 0.1 
Earth-like planets per Sun-like star. Habitable zones are expected 


to exist only around 0.6-1.4 M© {K to F5-class) stars (Kasting 
et al. 1993; Kopparapu et al. 2014), which make up 14.8% of stars 
by number (Chabrier 2003), so we take he = 0.015. 

Behroozi et al. (2013e) determined SFR(M*,t) for galaxies 
up to ~ 13 Gyr ago (z = 8), covering > 90% of all star formation 
(Behroozi et al. 2013c). The methodology is detailed in Appendix 
B; briefly, it involves linking galaxies at one redshift to galaxies 
with the same cumulative number densities at another redshift to 
trace their stellar mass buildup, as the most massive galaxies at one 
redshift will tend to remain the most massive galaxies at later red- 
shifts. The full computation also involves corrections for scatter in 
galaxy growth histories and galaxy-galaxy mergers (Behroozi et al. 
2013f,d). Knowing the stellar mass history of galaxies, one may use 
observed metallicity-stellar mass-redshift relations (e.g., Maiolino 
et al. 2008; Moustakas et al. 2011) or metallicity-stellar mass- 
star formation rate relations (e.g., Mannucci et al. 2010) to deter¬ 
mine galaxy metallicity histories (see also Munoz & Peeples 2015). 
Here, we use the fitting function in Maiolino et al. (2008), which is 
constrained for z < 3.5 (<11.7 Gyr ago), mildly extrapolated over 
the same redshift range as our star formation histories. As Maiolino 
et al. (2008) measure oxygen abundance ratios ([0/H]), we use the 


MNRAS 000, 000-000 (2015) 









On The History and Future of Cosmic Planet Formation 3 


Redshift 



Time Since Big Bang [GyrJ 


Redshift 



W) 

c 

’> 

o 


>> 


Q 


D 

C 

& 



Time Since Big Bang [Gyr] 


Redshift 

8 4 2 1 0.5 0.2 0 



Redshift 



Figure 2. Top-left panel: total Earth-like and giant planets formed in the Milky Way as a function of cosmic time. Giant planet counts have been shifted by 
a factor of 5 to allow better comparison with the Earth-like planet formation history. Top-right panel: average planet density in the Universe as a function 
of cosmic time. Earth-like planet formation tracks the galaxy/cosmic star formation rates, whereas giant planet formation times are greater at late times due 
to their metallicity dependence. Blue squares mark the median formation times of each population. The vertical dotted line indicates the formation time of 
the Solar System, which occuired after 80% of present-day Eai1h-like planets and 50% of present-day giant planets were formed in the Milky Way. Bottom 
panels: planet formation rates and densities, respectively, for the Milky Way and the Universe as a whole. Uncertainties in all estimates are I dex, arising 
from uncertainties in planet detection rates with Kepler. 


formula [Fe/H]= ^0.1 + 1.182 [O/H] (from fitting Milky Way stel¬ 
lar abundances in the SDSS-III APOGEE; Holtzman et al. 2015) 
to convert to iron abundance ratios. If all stars instead had solar 
iron-to-oxygen ratios, our derived giant planet abundance would 
increase by ~ 0.3 dex. 


3 RESULTS 

The resulting planet formation rates from Eq. I are shown in the top 
panels of Eig. I. Similar to star formation (Behroozi et al. 20l3e), 
planet formation rates per galaxy are greatest at early times in mas¬ 
sive galaxies. Indeed, Earth-like planet formation rates are exactly 
proportional to galaxy star formation rates (scaled by «£ = 0.0073) 
in our assumed model. The planet formation history of the Milky 
Way (present-day M* = 5 — 7 x lO*® Mq, from current Bayesian 
models; Licquia & Newman 2015) can be inferred by integrating 
PER(M*,f) along the median growth history for Milky Way-sized 
galaxies (dashed line, top panels of Fig. I). The planet formation 
history of the Universe, here expressed as the average volume den¬ 
sity of planet formation, is obtained by multiplying PFR(M*,f) by 


the volume density of galaxies as a function of stellar mass and 
cosmic time (0(M*,r), discussed in Appendix B). 

The product of PFR with (p shows when and where all Earth¬ 
like and giant planets formed (Eig. I, bottom panels). Typical 
galaxy masses at the time of planet formation are ~ M©, re¬ 

gardless of planet type. However, because of their metallicity bias, 
giant planets form later than Earth-like planets. Giant planet for¬ 
mation is rare in galaxies below lO^ M©; while the Magellanic 
Clouds may have some giant planets, it is unlikely that lower-mass 
dwarf satellite galaxies of the Milky Way will have any. In both 
cases, the lohnson & Li (20I2) minimum metallicity threshold is a 
weak one, as the vast majority of star formation has taken place at 
[Fe/H] > — 1.5. We find that total planet densities would be lowered 
by < 10% for Earth-like planets and ^ 0.01% for giant planets 
with this metallicity threshold (Eig. 1, bottom panels). However, 
this threshold would strongly diminish the number of Earth-like 
planets formed around stars older than 12 billion years in the Milky 
Way (Fig. 1), which is so far consistent with the age of the oldest 
observed star with Earth-like companions (11.2± 1.0 Gyr; Cam- 
pante et al. 2015). 

We show total planet formation histories and rates from Eq. 
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1 for the Milky Way and the Universe in Fig. 2. In the Universe’s 
observable volume (10*^ Mpc^), these results would imply ~ 10^® 
Earth-like planets and a similar number of giant planets.* Errors 
are dominated by planet incidence rates («£ and hq), which are un¬ 
certain at the 0.5-1 dex level (Lissauer et al. 2014) due to different 
detection efficiency estimates. Smaller systematic errors (0.2-0.3 
dex) come from uncertainties in initial mass functions, stellar popu¬ 
lation modelling in galaxies, and variation in individual galaxy star 
formation histories and metallicities (Behroozi et al. 2010; Peeples 
etal. 2014). 


4 DISCUSSION 

We discuss the Solar System’s relative formation time (§4.1) and 
its relation to the expected number of future civilisations (§4.2). 

4.1 Formation Time of the Solar System 

Eig. 2 shows that the Earth formed later than ~ 80% of similar 
planets in both the Milky Way and the Universe, matching previ¬ 
ous findings (Livio 1999; Lineweaver 2001). Comparatively, the 
Solar System (including Jupiter) formed closer to the median for¬ 
mation time for giant planets. This is not evidence for or against 
giant planets being prerequisite for life as there is a strong ob¬ 
server bias (Fig. 3). When calculating the age of our own planet, 
we are really calculating the time tc that it took our own species and 
civilisation to evolve. If U were extremely long, many new planets 
would have formed later than our own planet but before intelligent 
life evolved—so we would have concluded that our planet formed 
early compared to most other planets. However, as U is shorter than 
the current doubling time, tj, for stellar mass in the Milky Way 
(tc = 4.6 Gyr and tj ~ 20 Gyr), fewer planets have had time to 
form while civilisation has developed. Hence, the “late” formation 
time of our own planet speaks more to the ratio of tc to Ij than to 
conditions for habitability. 

This observer bias can be removed if we calculate our forma¬ 
tion time relative to all the planets which will ever be formed. The 
Milky Way is expected to merge with Andromeda (M31) in ~4 
Gyr (Cox & Loeb 2008), forming a single object with total (dark 
matter and baryonic) mass 3.17x10*^ Mq (van der Marel et al. 
2012). Using fitting formulae in Behroozi et al. (2013a) for its con¬ 
tinued mass growth, we expect that its total mass will asymptote to 
3.9 X10*^ Mq. Haloes of these masses are expected to have approx¬ 
imately the cosmic baryon fraction (16.5%; Hinshaw et al. 2013) 
of their mass in gas and stars (Werk et al. 2014), which translates 
to 6.4x 10** Mq of baryonic matter within the eventual halo. The 
present-day combined stellar masses of the Milky Way (Licquia & 
Newman 2015) and M31 (Tamm et al. 2012) are ~ 1.8 x 10** M©; 
correcting for the ~ 30% of stellar mass lost in normal stellar evo¬ 
lution (Chabrier 2003), this leaves 3.9x 10** M© of gas in the halo 
available for future star formation. As the remaining gas eventually 
cools and forms stars (as is expected to occur over the next trillion 
years; Adams & Laughlin 1997; Tutukov et al. 2000; Nagamine & 
Loeb 2004), this implies that the Earth has actually formed earlier 
than ~ 61% of all planets that will ever form in the Milky Way- 
M31 group. 


Slow Civilisation Formation 

, . , , _ , ..Median Age of Planets Formed Before 

Habitable Planet Formed*-, ; Civilisation Develops Telescopes 



DP Civilisation develops slowly, finds that own 

S planet formed before most existing planets 

Rapid Civilisation Formation (e.g., Earth) 

Habitable Planet Formed- 

^ Civilisation develops quicker, finds that own 

^ formed after most existing planets 

DP '• • Median Age of Planets Formed Before 

Civilisation Develops Telescopes 

Time-► 

■ Ongoing Planet Formation 
□ Time to Evolve Civilisation with Telescopes 

Figure 3. The relative formation time of one’s own planet depends on the 
time it takes one’s civilisation to form. As shown above, planet formation 
continues while civilisations are developing. Many planets will form if a 
civilisation is slow to develop, so by the time it is able to calculate its 
own planet’s formation time relative to others (~ the time when it develops 
telescopes), it will find that its planet formed early. In contrast, a rapidly- 
developing civilisation (e.g., ours) reaches that stage earlier, giving the Uni¬ 
verse less time to make more planets; the civilisation will then find that its 
own planet formed late relative to most others. 



Figure 4. Probability for the total number of planets with civilisations in 
the Universe, given that the Eaith formed before 92% of similar planets 
expected to exist. If Earth is the nth planet since the Big Bang to have 
formed a civilisation, then the average (expected) total number of planets 
with civilisations scales as 12.5n. Even for the most conservative possible 
assumption (i.e., that Eaith was the P* planet formed that evolved an intel¬ 
ligent civilisation), it is unlikely that we will be the only civilisation that 
the Universe will ever have {black line). As the number of earlier planets 
with civilisations increases (red and blue lines), it becomes more and more 
likely that the Universe will have many more civilisations than currently 
exist. For comparison, if the Milky Way today contained another civilisa¬ 
tion, it is likely that Earth would be at least the ten billionth planet to host 
a civilisation in the observable universe, which would eventually contain at 
least a hundred billion civilisations. 


^ The number of observable planets in our past lightcone, as opposed to 
planets currently in the Hubble volume, is ~ 10^^ for both Earth-like and 
giant planets. 
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Repeating this calculation for the Universe as a whole, we 
note that only 8% of the currently available gas around galaxies 
(i.e., within dark matter haloes) had been converted into stars at 
the Earth’s formation time (Behroozi et al. 2013c). Even discount¬ 
ing any future gas accretion onto haloes, continued cooling of the 
existing gas would result in Earth having formed earlier than at 
least 92% of other similar planets. For giant planets, which are 
more frequent around more metal-rich stars, we note that galaxy 
metallicities rise with both increasing cosmic time and stellar mass 
(Maiolino et al. 2008), so that future galaxies’ star formation will 
always take place at higher metallicities than past galaxies’ star for¬ 
mation. As a result, Jupiter would also have formed earlier than at 
least ~ 90% of all past and future giant planets. 

As shown in Fig. 2, planet formation rates have declined sig¬ 
nificantly since z ~ 2 (for Earth-like planets) and z ~ 1 (for giant 
planets), primarily because of declines in the cosmic star forma¬ 
tion rate. If these declines continue, most of the additional planets 
formed in both the Universe and the Milky Way will be in the very 
far future (100 Gyr to 1 Tyr from now) compared to the current age 
of the Universe (~ 13.8 Gyr; Hinshaw et al. 2013). Hence, as the 
Universe’s accelerating expansion is rapidly reducing the number 
of observable galaxies (Loeb 2002; Nagamine & Loeb 2003), most 
future planets formed in other galaxies will not be visible from the 
Milky Way. 


4.2 Probability of Other Civilisations 

The Drake equation (Drake & Sobel 1992) for calculating the 
number of intelligent, communicative civilisations (hereafter, just 
“civilisations”) is famously uncertain, with estimates of the civilisa¬ 
tion incidence per habitable planet ranging from 10^^ (Sagan 1963) 
to arbitrarily small values (e.g., < 10^^®, combining pessimistic es¬ 
timates from Carter 1983; Ward & Brownlee 2000; Schermer 2002; 
Spiegel & Turner 2012). Combined with our estimates of the num¬ 
ber of Earth-like planets (§3) and the fact of our existence, this 
would result in 1 to 10*^ civilisations in the Universe and 1 to 10'^ 
in the Milky Way at the present time. 

The formation time of our planet (compared to all which will 
ever form) gives weak but independent constraints on the total num¬ 
ber of planets with civilisations which will ever exist. Intuitively, 
if we were the only civilisation the Universe will ever have, the 
Copemican principle suggests that it is unlikely for our planet to 
have formed so early relative to other similar ones.^ As an example, 
we can calculate an upper bound for the chance that the Universe 
will only ever have a single civilisation (corresponding to civili¬ 
sation incidences of < 10^^* per habitable planet). For this upper 
bound, we adopt the prior that Earth is the first planet with a civil¬ 
isation to have formed—any possibility that Earth is not the first is 
incompatible with there being only one civilisation. We find using 
Bayes’ rule that the probability of there being N civilisations total 
is then: 

P{N\f,E = \)o.Nf-^P{H) (2) 

where / = 0.92 is the fraction of planets which have yet to form, 
£ = 1 is the assumption that the Earth has the first civilisation, and 


^ Our formation time relative to other habitable planets in the Milky Way is 
unexceptional (§4.1), giving little information on the total number of civili¬ 
sations our galaxy will have. Additional systematics for how planet migra¬ 
tion/scattering affect civilisation formation are discussed in Appendix A. 


P(N) is the prior on the number of planets with civilisations. Be¬ 
cause Eq. 2 falls off exponentially for large civilisation numbers N, 
it is only necessary to know the prior P{N) over a modest dynamic 
range (1 Al ^ 1000). The orders-of-magnitude uncertainties on 
parameters in the Drake Equation suggest that the prior P{N) on the 
number of planets with civilisations should be a log-normal distri¬ 
bution with an exceptionally large width (> 20 dex). Locally, then, 
it is an excellent approximation to take P{N) as uniform in loga¬ 
rithmic space (i.e., P{N) = j^) for 1 ^ A < 1000. 

The resulting probability distribution for the total number of 
planets with civilisations is shown in Fig. 4. The large fraction of 
planet formation which has not yet taken place (/ = 0.92) implies 
at most an 8% chance of us being the only civilisation the Universe 
will ever have. More typically, the expected total number of planets 
with civilisations would be {N) = 12.5. 

As noted above. Earth being the first planet with a civilisation 
is a very conservative assumption. For example, if the Milky Way 
today had another planet with a civilisation (~ 10^® civilisations 
per habitable planet), then Earth would be at least the ten billionth 
planet with a civilisation in the observable Universe. Generalising 
the problem, we suppose that Earth is the £th planet with a civili¬ 
sation. By Poisson statistics, the larger E is, the better the relative 
constraints on N will be. This is exactly analogous to an exposure 
time calculation: the more photons E that arrive in an 8s exposure, 
the better one can predict the total number of photons N if the ex¬ 
posure were extended to a full 100s. For planets, the generalised 
probability distribution is given by 

P(A|/,£)oc (A-£+l) (3) 

which is a binomial distribution with a prefactor (A — £ -f 1). As 
shown in Fig. 4, the expected total number of planets with civil¬ 
isations scales as (A) = 12.5£, and the relative uncertainties on 
this total drop as ~ . Hence, the more planets with civilisations 

which have formed before the Earth, the more likely it is for the 
Universe to continue forming many more in its future. 


5 CONCLUSIONS 

Current constraints on galaxy and planet formation suggest: 

(i) The Milky Way contains ~ 10® Earth-like and ~ 10*® giant 
planets (§3). 

(ii) The Hubble Volume contains ~ 10^® Earth-like planets and 
a similar number of giant planets (§3). 

(iii) A metallicity threshold of [Fe/H]= —1.5 has very limited 
effects on total planet counts (§3). 

(iv) Earth-like and giant planets both formed primarily in 
10*®'5 Mq galaxies; however, giant planets are much rarer than 
Earth-like planets in low-mass galaxies (§3). 

(v) Giant planets have median ages ~ 2.5 Gyr younger than 
Earth-like planets (§3). 

(vi) The Solar System formed after 80% of existing Earth-like 
planets (in both the Universe and the Milky Way), after 50% of 
existing giant planets in the Milky Way, and after 70% of existing 
giant planets in the Universe (§4.1). 

(vii) Assuming that gas cooling and star formation continues, 
the Earth formed before 92% of similar planets that the Universe 
will form. This implies a < 8% chance that we are the only civili¬ 
sation the Universe will ever have (§4.2). 
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APPENDIX A: ADDITIONAL SYSTEMATICS 
AEEECTING PLANET EORMATION AND STABILITY 

Beyond galaxy star formation rates and metallicities, planet forma¬ 
tion rates could vary with the stellar initial mass function (IMF) 
and the nearby stellar density. For the stellar IMF, host star mass 
is known to affect both the incidence and survival time of proto¬ 
planetary disks (see Williams & Cieza 2011, for a review). If the 
stellar IMF varied with redshift (Dave 2008) or galaxy mass (Wei- 
dner & Kroupa 2006), this would change the distribution of host 
star masses, correspondingly changing the planet formation rate per 
unit stellar mass. The Milky Way appears consistent with a non¬ 
evolving IMF (Chabrier 2003); more massive galaxies may have 
formed relatively more stars below 1 M© (i.e., a Salpeter 1955 
IMF) in the past (Conroy & van Dokkum 2012; Sonnenfeld et al. 
2012), although this has been debated (Smith 2014; Smith et al. 
2015). 

For Earth-like planets around sun-like stars, a Salpeter (1955) 
IMF yields only 30% more mass in F5-K stars than a Chabrier 
(2003) IMF. For giant planets, the incidence rate scales with host 
star mass to at least the first power for stars less than 1 M© 
(Gumming et al. 2008; Kennedy & Kenyon 2008), meaning that 
a Salpeter (1955) IMF would yield at most 70% more giant planets 
than a Chabrier (2003) IMF. These offsets are well within the 0.5- 
1 dex uncertainties on the overall planet incidence (Lissauer et al. 
2014). 

Stellar density may impact protoplanetary disk survival either 
via radiation (from nearby O stars) or via interactions in dense 
star clusters (Williams & Cieza 2011). Both would cause a red- 
shift evolution in the planet formation rate per unit stellar mass, 
because galaxies at high redshifts were (on average) much denser 
than galaxies today (van der Wei et al. 2014; Shibuya et al. 2015). 
For example, galaxies at z = 9 are ~ 20 times smaller in every direc¬ 


tion, resulting in galaxy densities similar to today’s globular clus¬ 
ters (Shibuya et al. 2015). Calculations for the effects of nearby 
O stars suggest, however, that the impact on protoplanetary disks 
is limited to regions beyond 15 AU (e.g., beyond Saturn) even at 
the dense centres of star clusters (Adams et al. 2004; Clarke 2007), 
with the effects limited to > 50 AU (e.g., beyond Neptune) for more 
typical environments within clusters (cf. Thompson 2013). For in¬ 
teractions, Wang et al. (2015) found that giant planet incidence in 
binary star systems was only reduced for binary separations less 
than 20 AU; incidence for binary separations between 20 and 200 
AU was instead mildly enhanced, and incidence at greater sepa¬ 
rations was similar to individual stars’ incidence rates. As a com¬ 
parison, stellar densities in globular clusters correspond to typical 
stellar separations of ~ 50000 AU; even for the longest-lived (10 
Myr) protoplanetary disks, cluster stars will have typical closest 
approaches of ~ 1000 AU (Adams et al. 2006). 

While planet formation rates may be relatively insensitive to 
the nearby environment, the same cannot be said of planet orbit sta¬ 
bility over several Gyr (see Davies et al. 2014, for a review). For¬ 
tuitous microlensing events have suggested significant numbers of 
free-floating giant planets (Sumi et al. 2011; Strigari et al. 2012), 
which could have scattered or migrated from their original birth¬ 
places but may also have formed in situ (Veras & Raymond 2012). 
In clusters, flybys of stars near planetary systems cause direct ejec¬ 
tion, longer-term (~100 Myr post-flyby) destabilization, and in¬ 
creased orbit eccentricities (Malmberg et al. 2011). That said, ob¬ 
servational evidence for different planet incidence rates in clusters 
has been mixed (Sigurdsson et al. 2003; Burke et al. 2006; Montalto 
et al. 2007; Quinn et al. 2012; Meibom et al. 2013). 

Accounting for these effects is beyond the scope of this pa¬ 
per. Instead, we note that typical stars in the Milky Way are 5-10 
Gyr old (Behroozi et al. 2013e), so that planets detected around 
these stars (e.g., with Kepler) are exactly those that have remained 
bound for long periods of time. Hence, the PFR is this paper is 
best interpreted as the bound planet formation rate. Even so, there 
is no guarantee that these planets remained in stable orbits for the 
lifetimes of their host stars; even mild changes in planet eccentric¬ 
ity could be detrimental to the development of civilisations (§4.2). 
Qualitatively, larger stellar densities in the past would lead to more 
interactions and therefore less planetary stability. This would re¬ 
duce the fraction of habitable planets formed before the Earth, pro¬ 
portionally raising the likelihood of future civilisations according 
to the argument presented in §4.2. 


APPENDIX B: RECOVERING GALAXY STAR 
EORMATION HI S TORIES 

Our reconstruction technique (detailed fully in Behroozi et al. 
2013e) uses forward modelling to extract the relationship between 
stellar mass, halo mass, and redshift (M*(M/,,z)). Briefly, we adopt 
a flexible parametrisation^ for M* (Af/,, z) ■ Any M* (M/,, z) in this pa¬ 
rameter space may be applied to a dark matter simulation, assigning 
galaxy stellar masses to every halo at every redshift. Linking haloes 
across redshifts with merger trees, the implied evolution of galaxy 
stellar mass, as well as average galaxy star formation rates, can be 

^ The z = 0 stellar mass—halo mass relation has six parameters: character¬ 
istic M* and M;,, a faint-end slope, a bright-end shape, a faint-bright tran¬ 
sition shape, and the scatter in M* at fixed M;,. For each z = 0 parameter, 
another variable controls the evolution to intermediate (z ~ 1) redshifts, and 
a third variable controls the evolution to high (z > 3) redshifts. 
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reconstructed. At the same time, the resulting predictions for ob¬ 
servables, including galaxy number densities, galaxy specific star 
formation rates, and total cosmic star formation rates are available 
from the assigned stellar masses and inferred star formation rates. 
Comparing these predictions to observations using a Markov Chain 
Monte Carlo method, we are able to constrain the allowable form 
of M*(M/,,z), and consequently the allowable reconstructions for 
galaxy star formation histories. 

Observational constraints are compiled in Behroozi et al. 
(2013e) from over 40 recent papers. At low redshifts, these include 
results from SDSS and PRIMUS (Moustakas et al. 2013); at high 
redshifts, these include constraints on galaxy number densities and 
star formation rates for z < 8 from Hubble observations (Bouwens 
et al. 2012, 2011; McLure et al. 2011; Bradley et al. 2012). Note 
that, above, we generate predictions for galaxy number densities 
as a function of stellar mass and redshift ((j){Mt,z))- For calcula¬ 
tions involving galaxy number densities in this paper, we use the 
prediction for (j){Mt,z) from the best-fitting M*(M/,,z); this gives a 
smooth form for (j){Mt,z) which is less susceptible to sample vari¬ 
ance than using the individual data sources in our compilation. 

We use the Bolshoi dark matter simulation (Klypin et al. 2011) 
for halo properties (including mass functions and merger rates). 
Bolshoi follows a periodic, comoving volume 250 Mpc on a 
side with 2048^ particles (~ 8 x 10®), each with mass 1.9 x 10* 
Mq, and was run with the ART code. (Kravtsov et al. 1997; 
Kravtsov & Klypin 1999) The adopted cosmology (flat ACDM; 
h = 0.7, = 0.27, CTg = 0.82, % = 0.95) is consistent with 

WMAP9 results (Hinshaw et al. 2013). Halo finding and merger 
tree construction used the ROCKSTAR (Behroozi et al. 2013b) and 
Consistent Trees (Behroozi et al. 2013d) codes, respectively. 
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